Intrinsic anomalous Hall effect in ferromagnetic metals 
studied by the multi-d-orbital tight-binding model 
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To elucidate the origin of anomalous Hall effect (AHE) in ferromagnetic transition metals, we 
study the intrinsic AHE based on a multi-orbital ((d xz , d yz )) tight-binding model. We find that a 
large anomalous velocity comes from the off-diagonal (inter-orbital) hopping. By this reason, the 
present model shows a large intrinsic anomalous Hall conductivity (AHC) which is compatible with 
typical experimental values in ferromagnets [10 2 ~ 10 3 fi - cm -1 ], without necessity to assume a 
special band structure at the Fermi level. In good metals where p is small, the intrinsic AHC is 
constant (dissipation-less) as found by Karplus and Luttinger. In bad metals, however, we find that 
the AHC is proportional to p~ 2 when h/2r is larger than the minimum band-splitting measured 
from the Fermi level, A. This crossover behavior of the intrinsic AHE, which was first derived in J. 
Phys. Soc. Jpn. 63 (1994) 2627, is recently observed in various ferromagnetic metals universally 



by A. Asamitsu et al. We also stress that the present (d x 
spin Hall effect in a paramagnetic state. 



, d yz )-tight binding model shows a huge 



PACS numbers: 72.10.-d, 72.80.Ga, 72.25.Ba 



I. INTRODUCTION 
A. Motivation and Purpose of the Study 

For a long time, the anomalous Hall effect (AHE) at- 
tracts increasing attraction from both theoretical and ex- 
perimental viewpoint. In general, the Hall resistivity is 
given by pn = R%B + AttR^M, where B is the mag- 
netic field and R° H is the ordinary Hall coefficient due 
to the Lorentz force. is the anomalous Hall coeffi- 
cient in the presence of magnetization M. In ferromag- 
nets, Rft is usually at least one order of magnitude larger 
than R? H . In paramagnetic heavy fermion (HF) systems, 
Ph/B takes a huge value due to the AHE since the uni- 
form susceptibility M/B — \ in HF is about 10 2 ~ 10 3 
times larger than that in usual metals owing to the strong 
Coulomb interaction. 

Recently, the AHE in transition metal ferromagnets 
has been intensively studied experimentally. For exam- 
ple, the AHE in ferromagnetic Pyrochlore oxides shows 
interesting behaviors [1-4]. The AHE is also observed 
in ferromagnetic spinels [5] and in Ru, Ti oxides [6, 7]. 
Several multilayer systems [8, 9] also show distinct AHE. 
The AHE in spin glass systems had been discussed in 
refs. [10, 11]. 

The theoretical study of AHE was initiated by Karplus 
and Luttinger (KL) [12] in 1954. They found that the 
anomalous Hall conductivity (AHC) a% y (= R a H M/ p 2 ) 
is finite and dissipation-less (i.e., a a y is independent of 
resistivity p) when M ^ 0. This KL-term is called the 
"intrinsic AHE" because it exists independent of impu- 
rities. In 1958, Smit presented a mechanism of "extrin- 
sic AHE" [13]: He found that spin polarized electrons 
are scattered asymmetrically around an impurity in the 
presence of spin-orbit coupling. The AHC due to this 
skew-scattering mechanism is linearly proportional to p if 



elastic scattering is dominant. The above two works were 
reproduced and refined in a unified way by Luttinger [14] 
by using the transport theory of Kohn and Luttinger, and 
his work were also reproduced by linear-response theory 
using diagrammatic technique [15, 16]. In 1970, Berger 
proposed another mechanism of extrinsic AHE, the side 
jump due to impurities [17]. This mechanism gives the 
AHC in proportion to p 2 . Note that the extrinsic AHE 
vanishes in a clean system without impurity. 

For a long time, the AHE had been generally regarded 
as an extrinsic effect, and the intrinsic AHE proposed by 
KL had been underestimated. This would be because the 
theoretical model assumed in ref. [14] was too oversim- 
plified. Moreover, ref. [14] could not offer a specific ex- 
pression for cr" Recently, however, various experiments 
suggest that the intrinsic AHE a xy oc (p)° is dominant in 
many transition metal ferromagnets [5, 18] even in good 
metals with p ~ lpQcm: It is surprising because the 
extrinsic AHC (oc p^ 1 ) should be dominant when p is 
small, if inelastic scattering is negligibly small. We here 
remind readers that the extrinsic AHC is very sensitive 
to the statistical property of impurity potential V; For 
instance, it vanishes when the average value of V 3 is zero 
[19]. 

After Luttinger's work [14], theory of intrinsic AHE 
has been developed by many authors. In general, a xy is 
composed of the "Fermi surface term" and the "Fermi sea 
term" ; the latter comes from quasiparticles inside of the 
Fermi sea, and it could exist even in insulators [20] . (Note 
that the conductivity <r xx is composed of only the Fermi 
surface term.) The term of the AHC derived by KL is a 
part of the "Fermi sea term" . Recently, M. Onoda and 
Nagaosa [24] and Sundaram and Niu [25] found that KL's 
AHC is expressed in terms of the "Berry curvature" [24, 
25] . (Note that Luttinger showed in his model that there 
is a term which almost cancels the Berry curvature term 
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[14].) On the other hand, Kontani and Yamada derived 
the intrinsic AHE due to the "Fermi surface term" based 
on the linear-response theory for the first time [26, 27]. 
They studied AHC in an orbitally degenerate periodic 
Anderson model, and succeeded in explaining the AHE 
in HF systems; o~ xy oc \ below the coherent temperature 
T , whereas o% y oc XP~ 2 above T [28, 29]. This results 
explain the experimental fact that the Hall coefficients in 
HF systems are proportional to p 2 below T . Miyazawa 
et al. proved that the intrinsic AHE also occurs in d- 
p models with orbital degree of freedoms [30]. Later, 
AHC's for Fe and SrRu03 were calculated based on the 
LDA band calculations [31, 32]. 

We still have to deepen the understanding of the mech- 
anism of AHE to explain experimental results in transi- 
tion ferromagnetic metals. For example, a recent exper- 
iment by Asamitsu et al. has revealed that the intrinsic 
AHC (a xy ~ 10 3 Q _1 cm _1 ) is observed in many ferro- 
magnets for p = 1 <~ lOO^ficm, whereas a xy starts to 
decrease in proportion to p~ n and n <~ 2 in bad metals 
where p S> 100/xficm. This drastic crossover is reminis- 
cent of the AHE in heavy fermion systems discussed in 
ref. [26]. Therefore, a detailed study based on an appro- 
priate model for transition metals is highly required. 

In the present paper, we study the intrinsic AHC in a 
tight-binding model with (d xz , d yz )-orbitals (e^-orbitals), 
which originate from the i2g-orbitals (d xz , d yz , d xy ) in the 
tetragonal crystalline field. The band structure of this 
model corresponds to a and (3 bands of Sr2Ru04 [30, 33- 
35] . We derive a general expression for the AHC valid for 
any damping rate h/2r, which enables us to study the 
AHC in bad metals. We show that the AHC is mainly 
given by the Fermi surface term for a wide range of h/2r 
since the KL's term is canceled by another Fermi sea term 
in a metallic state. We also find that the anomalous ve- 
locity due to d(xz)-d(yz) hopping gives rise to a large 
AHC comparable to experimental values (about e 2 /ha; 
a being the lattice constant), which will be a main ori- 
gin of a huge AHC in transition metal ferromagnets. In 
good metals, the intrinsic AHC is independent of r as is 
well known. However, it becomes proportional to p~ 2 in 
bad metals where H/2r is larger than the band-splitting 
around the Fermi level, A. This crossover behavior of 
the AHC can explain a recent observation of the AHE by 
A. Asamitsu et al. 



B. Origin of Intrinsic and Extrinsic AHE 

Here, we shortly explain the mechanisms of both the 
intrinsic AHE and the extrinsic AHE in more detail: The 
intrinsic AHE (both Fermi surface and Fermi sea terms) 
originates from the interband transition of quasiparticles 
due to off-diagonal terms of the velocity and the or- 
bital angular momentum I. (The diagrammatic expres- 
sion for o% y is given in fig. 4.) The off-diagonal velocity 
contains the "anomalous velocity" (e.g. v% oc k y ), which 
gives rise to the AHE. Since it is a purely quantum effect, 



it is difficult to find a simple classical analogue except for 
Rashba-type 2D electron gas model [36]. This might be 
a reason why the intrinsic AHE had been discounted for 
a long time. In the present model, a large anomalous 
velocity naturally comes from the inter-orbital (d xz -d yz ) 
hopping. By this reason, atomic d-orbitals are necessary 
to reproduce a huge AHC in real ferromagnets. 

On the other hand, the extrinsic AHE (the skew- 
scattering mechanism and the side-jump mechanism) 
happens even in a single-band model, although it is 
enhanced due to a multiband effect [13, 17]. Extrin- 
sic AHCs due to skew-scattering mechanism and the 
side-jump mechanism are proportional to p~ 1 and (p) , 
respectively. The former is caused by spin-dependent 
asymmetric scattering around the impurity in the pres- 
ence of spin-orbit coupling, and the latter is caused by 
lateral displacement of the wave-packet during the scat- 
tering [13, 17, 19, 37]. Diagrammatically, both skew- 
scattering and side-jump are expressed by current vertex 
corrections due to impurity potential, as shown in Ref. 
[19]. In general, detailed knowledge on the impurity po- 
tential is needed for a quantitative study of the extrinsic 
AHE. In the present paper, we do not study the extrin- 
sic AHE since recent several experiments suggest that the 
skew scattering term is small even in good metals [5, 18]. 

For a long time, it was believed that the AHE is in- 
dependent of the self-energy correction due to impurity 
potential. However, we show that the intrinsic AHC is 
reduced by the self-energy correction in bad metal. By 
this reason, change of the scaling law in the AHC at 
p ~ 100/if2cm reported by Asamitsu el at. [18]. is well 
explained in terms of the intrinsic AHE. 

II. MODEL AND HAMILTONIAN 

In the present paper, we study a square lattice tight- 
binding model with d xz and d yz orbitals, which is a sim- 
plified version of the (d xz ,d yz )-p z model. Miyazawa et 
al. [30] showed that the AHE exists in the latter model. 
However, they could not derive an explicit expression for 
AHC. Here, we derive explicit expressions for both the 
"Fermi surface term" and the "Fermi sea term" based on 
the present simplified model. 

The Zs-coupling term is indispensable for the AHE in 
ferromagnets: In a ferromagnetic metal with M || z, the 
Zs-coupling term Al • s is approximately proportional to 
—XMl Zl where A is the coupling constant. The AHE 
is caused by the inter-orbital transition of quasiparticles 
due to the off-diagonal elements of l z [12]. In a para- 
magnetic metal under B Zl on the other hand, Zeeman 
term for the angular momentum, ^B z l z , gives rise to 
the AHE [26]. Therefore, the AHE in paramagnetic met- 
als and that in ferromagnetic ones are caused by the same 
transport mechanism, although the origins of magnetiza- 
tions arc different. 

We consider that the present model with d xz - and d yz - 
orbitals describes a major part of the AHE in transi- 
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tion metal ferromagnets: In cubic or tetragonal crystals, 
(a\l z \/3) is nonzero only when (a, 0) = {xz,yz) [l z = ±1] 
and (a, (5) — (xy, x 2 — y 2 ) [l z — ±2]. Because levels of d xz 



and d yz orbitals are degenerate, the band structures com- 
posed of these orbitals are energetically close. This fact 
will be favorable to the intcrband hopping of quasiparti- 
cles between d xz and d yz orbitals, which is indispensable 

energy splitting be- 
^xy-orbital is of the order of 



J-yz 

for the AHE. On the other hand 
tween cL.2 _,.a -orbital and 



that a xy is proportional to (v x v v yy ), which can remain 
finite after the k-summation due to the anomalous veloc- 
ity. Therefore, a sizable AHC is caused by v x y and v y y . 
Atomic d-orbitals give rise to the huge AHC in transition 
metal ferromagnets. 

To realize the AHE in ferromagnets, the atomic Is- 
coupling is also necessary. It is given by H x = 



c k . Because \xz) = -{\l z = +l)-\l z = — l))/v / 2 



and \yz) = i(\l z = +1) + \l z = -1))/V2i, h> 



leV in square or tetragonal crystals. Therefore, the in- present basis is given by 



terband hopping of quasiparticles between d x 2 



and 



djy-orbitals would be difficult. As a result, the main con- 
tribution to the AHE in transition metal oxides will come 
from the inter-orbital transition between d xz - and d yz - 
orbitals. This is the reason why we study the (d xz , d yz )- 
orbital tight-binding model. 

Here, we represent the creation operator of an electron 
on xz (yz) orbital as (cj^). The Hamiltonian without 
Zs-coupling is given by H° = X^k^k^k^kj where 



hi 



Sk 

Sk Sk 



and 



t — ( r x\ „vt 
k — \ C \ 



= -2tcosk x , = 



(1) 



-2t cos k v 



and — At' sin k x sin k y . Here, —t and ±t' are the 
hopping integrals between nearest-neighbors and next- 
nearest- neighbors, respectively [34]. They are shown in 
fig. 1. 



d(yz) 
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\ 
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FIG. 1: (color online) (a) Hopping integrals between the same 
orbitals. (b) Hopping integrals between the different orbitals, 
which is given by the next nearest neighbor hopping. The 
sign of the hopping integral changes by 7r/2 rotation. This 
fact gives rise to the anomalous velocity. 



Then, the velocity matrix is given by 
where /i = x, y. They are given by 



2t sin k, r At' cos fc~ sin k„ 



At' cos k x sin k y 






At' sin k x cos k y 



dhl/dk^ 



, (2) 



At' sin k x cos k y 2t sin k v 



(3) 



We should stress that the off-diagonal elements of v x , 
v yx , is an odd- function of k y . In the same way, 



X 



,,yx 



is an odd-function of k x . They are called 
"anomalous velocity" . In later sections, we will show 



h x = sgn(s z )Af !/ , 



the 



(4) 



where A is the coupling constant and t y is the Pauli ma- 
trix for the orbital space. Note that l y = l z = in the 
present basis. In Fe, A = 70meV [31]. Hereafter, we put 
fi-Q = 1 for the simplicity of calculation. 



ocB 

G = 



a 



P 



a(3 
V v = 



...Jf 



FIG. 2: Diagrammatic expression for the Green function G a p 
and the velocity v" 13 . 

The Green function in the presence of atomic Is- 
coupling is given by Gk (u>) — (u> + [i — hf^ — h x ) ~~ 1 , which 
is given by 

G X ,X Gr 



xx VJ xy 

G yX Gyy 



d{uj) 



a k 



(5) 



where ak = £ k y + iAsgn(s 2 ) and d(uj) = (ui + /j, — £k)( w + 
k) - '~ ~ 



// — — |ak| 2 , which is expressed as 



E, 



k h 



(6) 



S + ^±V(^-^) 2 +%k| 2 , (7) 



{uj + fx — E^)(u) + fj, 
1 



d(oj) 



Et 



where £" k represents the quasiparticle dispersion. Fig- 
ure 3 shows the Fermi surfaces for (t,t') = (1,0.1). In 
Si'2Ru04, t'/t ~ 0.1 [33-35]. The electron density per 
spin, n, is set as 0.6 and 1.4. (Note that t,t' < in the 
present Hubbard model according to Slater-Koster [40]. 
However, we assume t, t' > because they are positive 
in Sr 2 Ru04 due to the presence of p 2 -orbital between 
the nearest i?it-sites [33-35].) The splitting A ± repre- 
sents the minimum band-splitting (\E^ — E^\) measured 
from the the Fermi surface of E^ -band. In Fig. 3, k* 
represents the position of the minimum band-splitting, 
A = min{A+, A"}. 

Here, we consider the quasiparticle damping rate 7 — 
K/2t, which is given by the imaginary part of the self- 
energy £k( e )- For simplicity, we assume that 7 is diago- 
nal with respect to orbital, and independent of momen- 
tum. Then, retarded and advanced Green functions with 
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FIG. 3: (color online) Fermi surfaces for (a) n — 0.6 and (a) 
n = 1.4. Here we put (t, t') = (1,0.1), which corresponds to 
Sr 2 Ru0 4 . 



finite 7 are given by 



G a p>(ijj - ij), 



(8) 
(9) 



where a,/? = xz,yz. This assumption will be justified 
when the damping is caused by local impurities since the 
local Green function g a ,p{c) — "■imp-f 2 X)k G Q ,^(k, e) is 
diagonal and independent of (a, (3) if A = 0. [Offdiagonal 
term of g a ,p(e) is proportional to A.] In the Born approx- 
imation, S(e) = rii m p/ 2 <7(e), where n; mp is the impurity 
concentration and I is the impurity potential. The Born 
approximation is valid when / is much smaller than the 
bandwidth IFband- In Appendix A, we derive the AHC in 
the Born approximation by taking the offdiagonal term 
of S(e) into account correctly, and find that its effect is 
negligible. 

On the other hand, the assumption given in eqs. (8) 
and (9) may not be valid in strongly correlated systems. 
In Appendix B and C, we study the intrinsic AHC when 7 
depends on momentum and band index, and find that the 
AHC could be changed prominently by these influences. 



III. ANOMALOUS HALL CONDUCTIVITY 

In this section, we derive the intrinsic AHC a" y (here- 
after, we drop the superscript a) and the longitudinal 
conductivity a xx for finite quasiparticle damping rate 7 
based on the linear-response theory [41]. Here, we take 
the self-energy correction due to correlation, which causes 
finite quasiparticle damping rate. On the other hand, 
we drop all the current vertex correction (CVC). The 
CVC due to impurities with Zs-coupling give rise to the 
skew scattering [19]. However, CVC due to local impuri- 
ties vanishes identically. Note that the CVC due to the 
Coulomb interaction does not cause the skew scattering 
[26]. Hereafter, we put the renormalization factor z — 1 
because z exactly cancels in the final formula of the AHC. 

From now on, we assume a complete ferromagnetic 
state with |-spin electrons only, i.e., sgn(s z ) = — 1. Here- 
after, we drop the factor e 2 /h = 2ire 2 /h (h being the 
Plank constant) in o^ t „ to simplify expressions. Accord- 
ing to Streda [20] , er M „ is given by 



' fay 



~I 1 „H 



(10) 



E 

k,aa'f3/3' 



de 



2ir " 



de 



/~tR r^R 1 /-1 A f~<A \ 



II 



\ E /^< Q ^/M 

k,aa'PP' 

® r<R n R r R ® n R 

w /~iA /~iA , i~iA w r^iA 

~~de ^ ' l3a ' /3 ' 9e l3a ' 



(11) 



(12) 



where a and G a p> are given by eqs. (2), (3), (8) and 
(9). /(e) is the Fermi distribution function. The dia- 
grammatic expression is given in fig. 4 (a). We call <7^„ 
and aj/y the "Fermi surface term" and the "Fermi sea 
term" , respectively, according to literature. 

The original KL's paper [12] studied only a part of the 
Fermi sea term (cr X y given in eq. (19)), which is equiv- 
alent to the Berry curvature term. On the other hand, 
Kontani and Yamada [26] derived the Fermi surface term 
a\. y in the J = 5/2 periodic Anderson model. They also 
showed that the Fermi sea term is absent in this model. 
In general models, however, both contributions exist as 
recognized by Luttinger [14]. Recently, a detailed anal- 
ysis for both terms was given in ref. [21] for graphene. 
However, comparison of each term was not complete for 
d-electron systems. In the present paper, we show that 
<j xy is much larger than cr* y in metallic systems for a wide 
range of 7 = H/2t. 

We note that the normal Hall conductivity due to the 
Lorentz force is proportional to 7~ 2 . In this case, a xy 
is much larger than cr x l ~ 0(7°) in good metals where 
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t = 1/27 1 is very large. On the other hand, the intrinsic 
AHC is ol the order ol 7 . By this reason, both <j\. y and 
all could be the same order. 



Xl 2 -M z 









(a) 



FIG. 4: (color online) (a) Diagrammatic expression for a xy . 
(b) First order term with respect to the atomic Zs-coupling in 
terms of the band-diagonal representation. 



From now on, we calculate the intrinsic AHC for 
(/i, v) = (x,y). First, we take the summation 
over a',a,/3',/3 in eqs. (11) and (12). After tak- 
ing the k-summation, only the terms (a', a, (3', (3) — 
(x,x,x,y),(x,x,y,x),(x,y,x,x),(y,x,x,x) remain fi- 
nite. Using the square lattice symmetry of the present 
model, we obtain that 



de J d R d A ' 



k v 
1, J 



(13) 



1 



1 

d^ 



(14) 



where v% x v* y 



8tt' sin 2 k x cos k y . Here, it is essen- 
tial that v% x Vy V is the totally symmetric representation 
(Ai g ), and it does not vanish after the k-summation. 
Note that the term proportional to G R G R + G A G A in 
a xy vanishes identically. 



The final result for the intrinsic AHC is given by 



'xy 



■xy 



Ila 
'xy 



lib 



I , Ila 

xy ' xy 



jib 

xy ' 



\ 

7T * — ' a 
k 



(( M -^) 2 +7 2 )((^-^ k ) 2 +7 2 )' 



2A 



1 



xlm 



k \~k ^k 

2/i - E+ - E k + 2» 7 
(A*-£^ + »7)(/*--Et+«7) 



1 



1^ n, xx v x V - 



x Im < In 



gk ~ A* - »7 



(16) 



(17) 



(18) 



(19) 



Below, we perform the numerical calculation for the AHC 
using eqs. (17), (18) and (19). 

In the same way, we calculate the longitudinal conduc- 
tivity for fi = v = x. In this case, (12) vanishes iden- 
tically, and only the Fermi surface term (11) remains fi- 
nite. After taking the k-summation and using the square 
lattice symmetry of the present model, we obtain the ex- 
pression for a xx at T = as follows: 



@X X 

la 

O XX 



a Ia _|_ a Ib 



(20) 



Jb 



^E{« 2; ) 2 ^-^ 2 +7 2 ] +2(^) 2 

k 

x[(r i ') 2 -A 2 + (/i-4)(M-^)+7 2 ] 
+4vTv*H xy (» ~ £„)} /d R (0)d A (0), (21) 

^Rc]T {{v%*f{fx £„ + z 7 ) 2 + 2(^) 2 

k 

x [(r 2 ') 2 - A 2 + ( M - ^ + 17) (A* - + n)] 

+4<x y r^ - e y + *7)> /(^(o)) 2 . (22) 



Here, we derive the expressions for a xy and cr^ at T = 
0. Then, the w-integration in eq. (14) can be performed 
using the relation 



dx 



(x~a) 2 (x-b) 2 



-(2[x-a-b) 
2 , I a — fx 



{a -by 



In 



/' 



(15) 



Note that a xx vanishes identically. In the case of 7 — > 0, 
a x % diverges as 7 _1 whereas a£x is finite. Therefore, 
&xx ~ m good metals. However, a xx and cr^| become 
the same order when 7 is very large. We perform the 
numerical calculation for the a xx using eqs. (21) (22) in 
later sections. 

Before performing the numerical study in §IV, we an- 
alyze eqs. (17), (18) and (19) when 7 is very small in 



6 



detail: In this case, 
1 



d R (0)d A (Q) 7 (E+ 
2^-E+-E k - 



Irn 



-^ k ) 2 +7 2 
2i7 



Im < In 



(p-E£+i<y)(ji-EZ 



(23) 



(24) 



- M - *7 

W-TT^-^+T^Oi-^). (25) 

Substituting above equations into (17)-(19), we obtain 
the following relation for 7^0: 



Ila 

'xy 

J xy 



2Xj2v x x x v x y y 

k 

S^-E+) + S(fi-E k ) 



'xyi 



k 



(26) 
(27) 

(28) 



According to eq. (28) 
comes from an area near k 
E k < are satisfied. When du '- 
(-e(»-E+) + B(»-E^))/d k 
In this case, 



the main contribution to al Ib 

y 



in fig 3, if E+ > and 
Ef£ — E k is small, then 



rib 



xx xy 



5(v-(E++E-)/2) 



(29) 



Note that this relation will not be satisfied in a spe- 
cial case where the Fermi level lies between a narrow 
bandgap. In the next section, we see that eq. (29) holds 
very well in the present model. 

Here, we summarize the obtained results in this sec- 
tion. 

(i) a J Jy was first recognized by Karplus and Luttinger 
[12]. It can be rewritten in terms of the summation of 
the Berry curvature [24] . 

(ii) en^ a is another "Fermi sea term", although its final 
expression (27) is written like a Fermi surface term. We 
find that aj y is very small in usual metals because a xy a 
and a xy b almost cancels each other. 

(iii) The Fermi surface term is canceled by a xy a in the 
clean limit, which is also recognized in 2D Dirac model 
[16]. However, this cancellation is imperfect when 7 is 
finite, as shown in the next section. 

Several previous works assumed that the intrinsic AHC 
is given by a x Ib [24, 31, 32]. One might consider that 



7^0. However, this assumption is not always guar- 
anteed as we will see in the next section: When 7 is as 
large as A, then (a) a xy w a xy still holds whereas (b) a xy 
becomes quite different from a xy . In Appendix B and 
C, we study the intrinsic AHC in the general case where 
7 is k-dependent or it depends on band index. In these 
cases, statement (a) and (b) are also true. Therefore, 
overall behavior of a xy for a wide range of 7 is well ex- 
pressed by the Fermi surface term a xy , not by cr xy b ■ This 
is an important result of this paper, which is obtained 
by considering all the terms contributing to the intrinsic 
AHC. 

Finally, we comment that both t' and h x change their 
signs under the gauge transformation \xz) — * — \xz). The 
AHC is invariant under this gauge transformation. 



IV. NUMERICAL STUDY 

In this section, we perform the numerical calculation 
for both <j xy and a xx at T = 0, assuming a complete 
ferromagnetic state where rij. — n and rij = 0. In this 
case, m z = ji^n. (We put /xb = 1 hereafter.) The main 
purpose of this section is to elucidate both the filling 
(n) and damping rate (7) dependences of the AHC. We 
perform the k-summations in eq. (16) for a xy and eq. 
(20) for <j xx numerically, dividing the Brillouin zone into 
5000x5000 meshes. 

The unit of conductivity in this section is e 2 /ha, where 
h is the Plank constant and a is the unit cell length (inter- 
layer distance in 2D systems). If we assume the length 
of unit cell a is 4A, e 2 /ha « 10 3 f2 _1 cm -1 . 




FIG. 5: (color online) Obtained total AHC (a, 
a function of n = m z . It is shown that a xy 



xy j and a xy as 
~ a xy is well 



satisfied, e /ha ~ 10 a f2 -1 cm -1 if we put the unit cell length 
a = 4A. o X y is not a monotonic function of m z , and changes 
its sign at m z ~ 1 

Figure 5 shows the n-dependence of a xy for n — 0.2 ~ 
1 v n = 2 corresponds to a ferromagnetic band insula- 



this assumption is justified by the present calculation for tors. In the present model, a xy (n) — —cr xy (2 — n). Here, 



7 



we put 7 = 0.05, which is sufficiently smaller than the 
minimum band-splitting measured from the Fermi sur- 
face, A, which is shown in Fig. 3. The bandwidth of 
the present model is approximately 4\t\ = 4. We put 
the is-coupling constant A as 0.1 ~ 0.2. If we assume 
t = 4000K, it corresponds to A = 400 ~ 800K, which 
are realistic values in transition metals. In the present 
study, \cr xy \ exceeds 0.5 x 10 3 fi _1 cm _1 at n ~ 1.8. The 
obtained magnitude of \cr xy \ is comparable with exper- 
imental value in Fe [31] and in SrRu03 [32]. Figure 6 
shows the A-dependence of \a xy \ for n — 0.6 and 1.4. 
\<7 xy \ is approximately proportional to A below A ~ 0.2, 
whereas it tends to saturate when A is as large as Aa=o 



since A 



E k , increases with A. [31]. 



In Fig. 5, a xy is positive (negative) for n < 1 (n > 1) 



a xy — at n 



1. The reason for the sign change of 
xy — is the following: According to eq. (16), the main 

comes from an area 
tt' sin 2 fc* cos k* . 



a 1 or a IIb ) 

xy w xy I 



contribution for a xy 
near k* . Therefore, a xy ~ v x Xv y V \ v , 
Because tt' > 0, the sign of a xy changes to be negative 
when |fe*| = |fc*| exceeds tt/2. The sign change of the 
AHC as a function of m z is actually observed in SrRuC>3. 



t'/t=0.1, y=0.05 




tional to 7 2 for \E* — p\ -C 7 in J = 5/2 periodic Ander- 
son model (PAM). (E* is the energy of /-electrons.) The 
different 7-dependence of a xy between in d-p model and 
in PAM comes from the fact that the conduction band- 
width Wband ' s mucn wider than the heavy quasiparticle 
band in PAM. In the PAM, the quasiparticle damping is 
given by 7 = T ■ V 2 /((p, - E^f + T 2 ), where V is the c-f 
mixing term and T is the imaginary part of the /-electron 
self-energy. Because 7 saturates at T ~ \fj, — E^\ and the 
relation V 2 /{n - E?) ~ W^ and holds in the PAM, 7 will 
always be smaller than W£ and even if T ^> W^ and . 
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FIG. 7: (color online) 7-dependences of a xy (= cr tot ), al y (= 
a 1 ), ai y a (— a IIa ) and ai y b (= a IIb ), respectively. a xy shows 
a crossover behavior at 7 ~ A. It is shown that a^ y repro- 
duces a overall behavior of the total a xy for a wide range of 
7- 



FIG. 6: (color online) A-dependence of a xy for n = 1.4 and 
0.6. 

Next, we examine the 7-dependence of the AHC. Fig- 
ure 7 shows the total AHC (cr xy ), a xy , —c^iy " and a 1 ^ for 
n = 1.4. We see that all of them are 7-independent when 
7 -C A. On the other hand, they start to decrease with 
7 when 7> A: This is easily recognized from the func- 
tional forms of eqs. (17)-(19). We see that a xy is almost 
equal to —&i y a and a xy w a 1 ^ when 7 <§C A, as discussed 
in the previous section. On the other hand, a xy w cr l x%j 
whereas a xy is quite different from a X y b for 7 ^> A. As a 
result, the Fermi surface term a xy succeeds in reproduc- 
ing the overall behavior of the AHC for a wide range of 
7. <T xy is proportional to 7 -3 for 7 ~ Wband, which can 
be recognized by eq. (17). 

This crossover behavior of a xy at 7 ~ A was first 
pointed out by ref. [26]. They found that a xy is propor- 



We also study the 7-dependences of the longitudinal 
conductivity <r xx . Figure 8 shows that the relations 
(Jxx oc 7 _1 and a x % 3> <r x b x hold well for a wide range of 7. 
The conductivity a xx does not show any clear crossover 
behavior at 7 ~ A. When 7 ~ IFband, o'xx seems propor- 
tional to 7~ 1,5 , and cr^ is as large as a xx . Although eqs. 
(21) and (22) suggest that a xx cx 7 -2 when 7 ^> Wband, 
we could not find such a behavior for 7 > Wband- 

Finally, we study the relation between the AHC and 
the resistivity p = l/a xx . Figure 9 shows that a xy is in- 
dependent of p when p < 0.1 [~ 100/ificm]. On the other 
hand, a xy starts to decrease for p > 0.1 in proportion 
to p~ 2 . This crossover behavior of a xy at p ~ 100/ificm 
is observed universally by recent experiments on various 
transition metal ferromagnets by Asamitsu et al.[18]. 
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FIG. 8: (color online) 7- dependences of <r^, a xx and a xx = 



' XX < ^ xx • 
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FIG. 9: (color online) Obtained relation between a xy and p = 
for n = 1.6 and 0.4. p — 0.1 corresponds to 100/^ficm. 
It reproduces well the experimental universal "Asamitsu plot" 
for various transition metal ferromagnets, where a xy = 10 2 ~ 
cm for p = 1 ~ 100/if2cm, whereas a xy oc p for 
p > lOOpficm. 



V. ANALYSIS OF ANOMALOUS HALL 
CONDUCTIVITY IN THE BAND-DIAGONAL 
REPRESENTATION 



In this section, we analyze the intrinsic AHC in terms 
of the band-diagonal representation, which has been ana- 
lyzed by many authors. We note that this representation 
is useful for the analysis of the AHC only when 7ct is 
orbital independent. In this representation, according to 



eqs. (11) and (12), 
1 , 



E 



-V' VI 



(30) 



TI 



^ lie 

" fc ml Im 

^ I 2?r x v 



xlm 



1 



(e-El + z 7 )2 (e - E£ + i 7 ) 



(e-El + i 7 ) (e 



EP 



(31) 



at zero temperature. Here, I and m are the band suffices, 
and we dropped the diagonal terms I = m because they 
vanish identically. After performing the e- integral in cr^, 
the Fermi sea term is expressed as a*.y a + a 1 ^: 



Ila 
'xy 



rib 



2lT " x " v E L~ E ™ 



xlm 



K + E k - 2 ^ - 2i l 
(4- M + -; 7 )(£-- M -i 7 ) 



,(32) 



v E < l <' 



x Im < In 



(ei - E™r 



E L - m - «7 



El; 



fi-vy 



(33) 



When 7 is very small, they are given by 



■nj 



Ila 

'xy 

lib 

xy 



E *0*-4)- 



ml Im _ Im ml 



K ~ E F 



-a 



xy 



* E /(^ 



ml Im _ Im ml 
I \ U x u y u x v y 



(34) 
(35) 

(36) 



jlm/^m _ ^ and jln 



Because v.™ 

(k, l\d/dk^\k, m), a 1 ^ is equivalent to eq. (2.16) of ref. 
[12], or eq. (3) of ref. [24]. Based on eq. (36), When 
the Fermi energy lies in a gap, a 1 ^ is quantized when 
the Fermi energy lies in a gap [42]. In the metallic state, 
on the other hand, a 1 ^ is rewritten as the Fermi surface 
term due to the partial integral. 

Off course, we can reproduce eqs. (26)-(28) by 
applying eqs. (34)-(36) to the present tight-binding 
model given in §11. In the band-diagonal representation 
(E£; a = ±), we can show that the off-diagonal velocity 
is given by 

v+- = -{\ a \ 2 vT-c y (r-e)v x x v 

iXv^\E+-E^\}/(\a\\E+-E^\), (37) 

«+- = {\a\\f-e v {e-i v K v 

-i\v?\E+ - E^\) /(\a\\E+ E£\), (38) 
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and + = {w+ }*. Therefore, 



xx xy i xy yy 
aa aa ■ \ x V x u y 

< Vy = -*A L _ F a 

+ (non Ai g terms), 



(39) 



where a = ± and a — —a. Substituting this result into 
eqs. (34)-(36), we obtain the same results given in §111. 



VI. DISCUSSIONS 
A. Summary of the Present Study 

In the present paper, we studied the mechanism of the 
intrinsic AHE in transition metal ferromagnets based on 
the (d xz , c?j, z )-orbital tight-binding model. The origin of 
the anomalous velocity is the hopping integral between 
d xz and d yz orbitals, t' . By virtue of it, we could re- 
produce a typical experimental value of the AHC in fer- 
romagnets; 10 2 <~ 10 3 O _1 cm _1 . Thus, the anomalous 
velocity due to atomic d-orbitals will be the main ori- 
gin of the AHE in transition metal oxides. This fact 
has been overlooked in previous theoretical works based 
on electron gas models without atomic orbitals. We 
note that the present model breaks the parity symme- 
try [(x, y) — > (— x, y)] due to t' ^ 0. This is a necessary 
condition for a spontaneous Hall effect [24, 44] . 

In accord with the present study, paramagnetic com- 
pound Cai.7Sro.3Ru04 shows large AHE under the mag- 
netic field [43]. Its magnitude is comparable with the 
large AHE in /-electron systems (such as UPt 3 ) which 
originates from the anomalous velocity due to atomic /- 
orbitals [26]. The Fermi surfaces of this compound are 



composed of ^-orbital (d xz , d yz , d^-orbitals); 7-band 
composed of d^-orbital is absent in the present {d XZl 



-tight-binding model. We will study the AHE of this 
compound in more detail based on the ^-orbital tight- 
binding model, which reproduces the bandstructure of 
(Ca,Sr)2RuC>4 accurately [23]. 

We derive a general expression for the AHC valid for 
any quasiparticle damping rate 7, by performing the an- 
alytic continuation carefully. [Equations (17), (18) and 
(19) for (d xz , dj, z )-orbital model, and eqs. (30), (32) and 
(33) for a general model.] Using the general expression, 
we succeeded in explaining the experimental crossover 
behavior of the AHC in bad metals; a xy is constant for 
p < 100/if2cm, whereas a xy oc p~ 2 for higher resistivity. 
This overall behavior is mainly given by the Fermi sur- 
face term, <J xy7 whose importance was intensively studied 
in ref. [26]. 

We stress that the intrinsic AHC in the present model 
is not a monotonic function with respect to m z = n j . In 
partial ferromagnets, the total AHC is given by 



cr 2;y (n T ,n i ) = cr xy (pi) - <J xy {^) 



(40) 



the sign-change of a xy occurs in correspondence with sign 
of v* v at k*. We expect that this result will explain the 
sign-change of the AHE in SrRu03 as a function of the 
magnetization. We comment that authors of ref. [24] 
found that a large value of AHC appears from a xy b when 
Fermi level lies inside a narrow " anticrossing band gap" , 
and they suggested that this mechanism accounts for a 
huge AHC in ferromagnets. However, this condition will 
be satisfied only for a narrow range of magnetization. 
The present model with (d xz , d yz )-orbit&\s can give an 
enough magnitude of the AHC, although band crossings 
are absent (except at k = (0,0) and (it, it)). We stress 
that the AHC in the present model shows a sign change 
with respect to m z , even if the topology of the Fermi 
surface is unchanged; see Fig. 5. 

The relation a xy — <r xy b (Berry curvature term) will 
not hold in usual metallic compounds: For example, in 
bad metals where 7 is as large as the minimum band- 
splitting measured from the Fermi surface, A, the rela- 
tions a xy i=a a xy (Fermi surface term) and \a xy \ ^> Wi y \ 
hold well for a wide range of 7 as shown in fig. 7. Even in 
good metals, the relation a xy = <r xy b is also invalid when 
(a) 7 is k-dependent, and/or (b) 7 depends on band in- 
dex, i.e., 7; /7m 7^ 1: In these cases, the Fermi surface 
term deviates from eq. (26) as shown in Appendix B and 
C. 

In summary, the relation a xy = <r xy b will hold only in 
good metals as well as when 7;(k, ui) is independent of I, 
k and u> On the other hand, the relation 

a xy ?y a xy (Fermi surface term) (41) 

will be universal in real metallic systems, since the Fermi 
sea terms almost cancel each other except for a special 
situation where p lies inside a narrow anticrossing band 
gap. This result is consistent with recent theoretical work 
on metallic graphene [21]. For a quantitative study of the 
intrinsic AHC, however, we have to calculate <r xy , a xy a 
and al y b on the same footing. This is also important 



where cr xy (p) is given by eq. (16), and Pm) is the chemi- 
cal potential for the electron density n-my We found that 



'xy 

conclusion in this paper. 

Finally, we comment that the present (d xz , d yz )-tight 
binding model in a paramagnetic state shows a finite 
spin Hall conductivity (SHC) a xy , that is, cr z -spin current 
along y-axis occurs under the electric field along x-axis. 
a z xy is given by (— ft/e) times eqs. (16). As shown in Fig. 
5, a z xy reaches 0.6 P/e)(e 2 /ha)] - 600 [he- 1 ■ Vl^cmT 1 ] 
for a =4A, which is almost one order of magnitude larger 
than typical SHC in semimetals [22]. In later publica- 
tions, we will study the spin Hall effect in more detail 
[23] 



B. Comparison with Experiments: Transition 
Metals 

Asamitsu et al. examined the experimental relation 
between the AHC and the resistivity p in various fer- 
romagnets. They observed the dissipation-less intrin- 
sic AHC in the intermediate conducting region with 
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p = 1 ~ 100/if2cm. In the bad metal region, a xy is pro- 
portional to p~ 2 , which is well explained in the present 
study. We found this crossover at ~ 100/i£lcm occurs 
when 7 ~ A. According to Fig. 8, lOO^ficm corresponds 
to A ~ 0.2 w 800K. Another crossover to the good metal 
region (<~ 1/zQcm) seems to be rather complex: It might 
be due to the extrinsic AHE or a phenomenon related to 
the normal Hall effect. 

Recently, S. Onoda ct al. [45] calculated both the in- 
trinsic AHC [12, 26] and the extrinsic AHC [13, 19] in 
a Rashba-type 2D electron gas (2DEG) model. They 
claimed that the extrinsic-intrinsic crossover occurs when 
7 <~ A, and tried to explain the crossover at p <~ 1/zQcm 
observed by Asamitsu. However, it will not be true be- 
cause 1/iOcm corresponds to 7 ~ 10~ 3 w 8K, which is too 
small for a realistic value of A. (see Fig. 8.) We stress 
that Inoue et al. [46] proved that the intrinsic AHE in 
a Rashba model almost vanishes due to the cancellation 
by the current vertex correction (CVC) by impurities. 
On the other hand, we have checked that the CVC for 
anomalous velocity is absent in the present (d xz ,d yz )- 
orbital model within the Born approximation, when the 
impurity potential is delta-functional. By this reason, the 
intrinsic AHC obtained in the present work is not mod- 
ified by the CVC. We conclude that the present model 
with atomic d-orbital degrees of freedom is essential for 
a realistic study of AHE in ferromagnets. 



large intrinsic AHC [26, 27]. They predict that a xy oc \ 
below T , whereas <j xy oc XP 2 above T . We expect 
this intrinsic AHE is widely observed since the J = 5/2 
PAM is a well-established effective model for Ce-based 
HF systems: For example, it could explain a large Van- 
Vleck susceptibility in Kondo insulator [47, 48] . 

There remain a lot of future works to be done about the 
AHE in HF systems. For example, Nakajima et al. found 
that the AHC in two-dimensional HF system, CeMIns 
(M=Co,Rh,Ir) is negligibly small [49, 50]. In this com- 
pound, strong temperature dependence of i?n is given by 
the normal Hall effect, due to the effect of the current 
vertex corrections (or backflow) in nearly antiferromag- 
netic Fermi liquids [51, 52]. In addition, an interesting 
field-direction dependence of AHC is found in CeCug or 
CeCus.gAuo.i [53]. Theoretical studies on these experi- 
mental results will give us important hints to understand 
the electronic states. 
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C. Comparison with Experiments: Heavy Fermions 

In paramagnetic heavy fermion (HF) systems, the Hall 
coefficients takes huge values due to the AHE, since 
the uniform susceptibility M/B = x in HF is promi- 
nently enhanced by the strong Coulomb interaction. In 
HF systems, the large quasiparticle damping rate h/r 
comes from the strong electron correlation (or Kondo 
resonance). In the early stage, Coleman et al [38] and 
Fert and Levy [39] developed theories of the extrinsic 
AHE: developed theories of extrinsic AHE: They studied 
the extrinsic mechanism based on the /-electron impurity 
Anderson models with rf-orbital channel, and predicted 
the relation i?H oc XP above T . They assumed that the 
impurity Anderson models is applicable for the study of 
AHE in periodic systems if T > To. 

On the other hand, Kontani and Yamada found that 
the J = 5/2 periodic Anderson model (PAM) gives a 



APPENDIX A: AHC WITHIN THE BORN 
APPROXIMATION; EFFECT OF OFFDIAGONAL 
ELEMENT OF 7 

In previous sections, we studied the AHC by assum- 
ing that ImS is diagonal with respect to atomic orbital. 
This assumption seems to be justified since llmEi^l » 
|ImSi.2| when Zs-coupling constant A is much smaller 
than bandwidth Wband- [When A = 0, ImEi^ = since 
the local Green function g(w) = Y] k Gk(u;) is diagonal 
due to the orthonormality of rf-orbitals.] Here, we de- 
rive the AHC within the Born approximation, by taking 
both the diagonal and off-diagonal elements of ImS. As 
a result, it is confirm that the off-diagonal element is 
negligible in calculating the AHC in the present model. 

In the Born approximation, the Green function is given 

by 



rR(,,\ - 1 ( u + p-£l + il al + + g - ih) , , 

° [U) ~ d{u) \a£-i\{l+g-ih) w + /i -g+i 7 ^ 

G» = -L ( £ + - ^ " J? < +; A (! +9 + ih ) ) , (A2) 

v ' d(u) \ <*%. ~ «A(1 + g + ih) u + /z - - 17 J' y > 



I 

where 7, g and h come from the self-energy correction in the self-inconsistent Born approximation: 

i-y = -TmS^(O) 

= -^n im p7 Im\_^^p 
i\(g-ih) = £f 2 (0) 
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where n; mp is the impurity concentration and / is the 
impurity potential. Therefore, 7 = 7m imp / 2 iV(0)/2 ~ 
ni m pl 2 /Wband, where N(0) is the DOS at the Fermi 
level. In the same way, we can estimate that g, h ~ 



(i/w hand y ~ 7 /wi 



band 



< 1. 



We can derive the expression for the AHC by inserting 
eqs. (Al) and (A2) into the Streda formula in §111, in 
the same way of deriving eqs. (13) and (14). Using the 
square lattice symmetry of the present model, the Fermi 
surface term at zero temperature is given by 



2A 



^ d R (0)d A (0) 



[-Y(l + g) + h( Ji -?)] 



(A5) 

Since the main contribution in the k-summation comes 
from k ~ k*, where k* represents the position of the 
minimum band-splitting A, as shown in fig. 3. Since 
tyl* ~ £ y \ ~ 7 at k = k*, the second term of eq. 

(A5) is negligible. Therefore, the correction due to off- 
diagonal element of the self-energy correction, i\(g — ih), 
is very small. As a result, eq. (13) is justified in the Born 
approximation. 



o 



o 





(b) 



FIG. 10: (a) Diagrammatic expression for the intrinsic AHC 
given in eq. (A5). Self-energy correction due to impurity 
scattering is not diagonal when A / 0. (b) and (c) represent 
the Born and the second Born approximation, respectively. 
The latter gives the skew scattering extrinsic AHE when the 
impurity potential is non-local, (d) represents the side jump 
extrinsic AHC. 

Figure 10 (a) represents the intrinsic AHC given in eq. 
(A5). In good metals, the diagonal self-energy correction 
does not affect the AHC, and the AHC due to off-diagonal 
self-energy correction is negligible as we discussed above. 
In addition, current vertex correction (CVC) term (b) 
also exists in the Born approximation. Figure (c) and 
(d) represent other CVC terms which gives the extrinsic 
AHC: (c) gives the skew scattering extrinsic AHC that 
is proportional to a xx , and (d) represents the side-jump 
extrinsic AHC term that is proportional to {cr xx }°, re- 
spectively. The latter is caused by the lateral displace- 
ment of the wave-packet during the scattering, which is 
theoretically expressed by the modification of velocity op- 
erator due to impurity potential [17, 19]. In the present 
model, (b)-(d) vanish identically in the case of local im- 
purity potential. In contrast, extrinsic AHC occurs by 
local impurity potential in 2D electron gas models when 
the Zs-coupling term is k-dependent [16]. 



APPENDIX B: AHC WHEN 7 IS MOMENTUM 
DEPENDENT 



In §111, §IV and §V, we have neglected the momen- 
tum dependences of the damping rate 7 to simplify the 
discussion. However, this assumption is not realistic be- 
cause lm£k(0) is usually k-dependent. In this appendix, 
we study the AHE by taking the k-dependence of 7 into 
account. Then, we have to distinguish j x — ImE xx (— iO) 
and 7 y = lmY, vv (— iO). We show that the expression for 
a^y in eq. (26) could be changed prominently when the 
k-dcpendence of 7 is taken into account. 

Along with the derivation of eq. (13), we obtain the 
expression for a xy in the present case: 



a ^ v% x v x y y lv + v*y v yy 7; 



d R (0)d A (0) 



(Bl) 



where 



d R (0)d A (0) = ((p E+f + 7 »)(0i - E k f + 7 »).(B2) 
When both j x and j y are very small, 



7 



2 V' E+-E^ lx + ly 



(B3) 



In this case, eq. (Bl) is given by 

I v x x Vy V Jy + v x x v v y y y lx 

S^-E+^+S^-E^) — 
7 k 7 k . 



(B4) 



When 7^ is momentum dependent, j x ( y )/~f + , lx( y )/l~ 7^ 
1 except for | fc^ | = \k y \. Therefore, eq. (B4) is not equal 
to eq. (26) even in good metals ('J x ,jy — > 0). 

In general, the momentum dependence (and the band- 
dependence) of 7 could give rise to a nontrivial modu- 
lation of the AHC due to the Fermi surface term, cr^„. 
This effect will be prominent in strongly correlated elec- 
tron systems like high-T c cuprates and heavy Fermion 
systems. On the other hand, the Fermi sea term cr^, 
which is given by eqs. (27) and (28), is satisfied when 
lx, ly -> 0, even if 7^/7;, 7^ 1. 



APPENDIX C: CALCULATION OF AHC WHEN 
7 IS BAND DEPENDENT 

As shown in Appendix B, we have shown that the ex- 
pression for Fermi surface term a xy given in eq. (26) is 
modified by the k-dependence of 7. In this section, we 
explain that eq. (26) is also changed when 7 depends on 
bands. This fact was already pointed out in Appendix 
D of ref. [26]. Here, we derive the first order term of 
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T 1 

'xy 



cr*„ with respect to the /s-coupling, and correct a mis- 



take in Appendix D of ref. [26]. By extracting a single 
^s-coupling term from eq. (30), we obtain 



■ xy 



2tt 



E G ? G ? 

x [v l ^G^vf(l\l z \m)+Giv l y m (m\l z \l)) 
+(v l x m G^(m\l z \l)+vg l G^(l\l z \rn))v l y l ] (CI) 



where Gf" = (/i — E l k + i'ji)^ 1 and v ™ are the Green 
function and the velocity without Zs-coupling term. The 
diagrammatic expression is given in fig. 4 (b). Using the 



relation (l\l z \m) 
that 



- E 

2tt ^ 



(m\l z \l) and v l ™ — v™ 1 , we obtain 
" G?Gt(G*-G*)(l\l x \™) 



When 7; is very small, 



(C2) 



2m 



(C3) 



In this case, eq. (C2) is given by 



^E 



k,Z^m 

i(l\l z \m) 
; (El - Egf 



[v l H m ~ v ^ v y] ■ ( C4 ) 



In summary, a xy depends on the ratio of j m /ji- &iy 
could take a huge value if ji differs much from band to 
band. This fact was already pointed out by ref. [26]. 
On the other hand, eqs. (27) and (28) for the Fermi sea 
term a^y i s satisfied when ji — > 0, even if l m /ll 7^ 1- 
The present result will be valid even if the Is-coupling is 
treated unperturbatively, as done in §11 - §IV. In fact, 
Fig. 6 shows that <j xy is approximately linear-in-A for 
A<0.2. 

In Appendix D of ref. [26], Eg in 8{p, - Eg) was 
replaced with E^ by mistake. Thus, eq. (C4) is not equal 
to the KL's term which is given by [eq. (3.5)]x(2/m^) 
in ref. [12] with (a, b) = (x, y). This difference is natural 
since the former and the latter are the Fermi surface term 
and the Fermi sea term, respectively. 
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